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We study transport of interacting particles in weakly disordered media. Our one-dimensional 
system includes (i) disorder: the hopping rate governing the movement of a particle between two 
neighboring lattice sites is inhomogeneous, and (ii) hard core interaction: the maximum occupancy 
at each site is one particle. We find that over a substantial regime, the root-mean-square displace- 
ment of a particle, a, grows super-diffusively with time t, a ^ (et)^^^, where e is the disorder 
strength. Without disorder the particle displacement is sub- diffusive, a ~ t^^'^, and therefore dis- 
order dramatically enhances particle mobility. We explain this effect using scaling arguments, and 
verify the theoretical predictions through numerical simulations. Also, the simulations show that 
disorder generally leads to stronger mobility. 
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Disorder, inhomogeneities, and impurities are ubiqui- 
tous in physical systems and are widely used to control 
properties of matter. Some of the most fascinating phe- 
nomena in contemporary physics including localization 
[Hi SSI J glassiness 4, 5], slow relaxation 1^6,], and frustra- 
tion [7[ are unique consequences of disorder. 

While the effects of disorder on noninteracting parti- 
cles are well-understood, the consequences of disorder on 
interacting, strongly correlated particles remain an open 
question (11,01 . In a quantum system, an isolated particle 
is localized by disorder, but localization can be destroyed 
when there are two interacting particles 13, Hi- Hence, 
disorder and particle interactions compete. We investi- 
gate this interplay between disorder and interaction in a 
classical system where inhomogeneities are known to trap 
particles and drastically decrease their mobility, and find 
that, as in the quantum case, disorder has opposite effects 
on noninteracting and on interacting particles. Our main 
result is that disorder speeds up the motion of interact- 
ing particles whereas disorder slows down the movement 
of noninteracting particles. 

We generalize the standard exclusion process 13, l3| 
to study the interplay between disorder and interactions 
II E H Ezl- Our system is an unbounded one- 
dimensional lattice whose sites may be either occupied 
by a single particle or vacant. Initially, the lattice is 
populated at random by identical particles with concen- 
tration c. Each particle may hop from an occupied site 
into a neighboring vacant site and this diffusion process 
is governed by the following rates: p+{i) is the hopping 
rate from site i to site i + I, and similarly, P-{i) is the 
hopping rate from site i to site i — 1. While the total 
hopping rate is uniform, and is set to one without loss 
of generality, p~{i) = 1, the lattice is inhomo- 

geneous. At every site there is, with equal probabilities, 
a bias to the right, p+ = 1/2 -I- e, or a bias to the left, 
p+ = 1/2 — e, as illustrated in figure [TJ The parameter 
0<e<l/2is the disorder strength. Note that the disor- 
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FIG. 1: Illustration of the disordered interacting particle sys- 
tem. The arrows indicate the bias at each site, the circles 
indicate vacant sites, and the bullets indicate occupied sites. 



der is quenched, uncorrelated, and uniform in strength. 
Moreover, since every lattice site accommodates a sin- 
gle particle, the particles interact via hard core repulsion 
[l8|. Our problem gen eralizes two well-known processes: 
single-file diffusion ^1 [111, [H, [H, [H [H [ll with in- 
teraction but no disorder corresponds to the limit e ^ 
, and Sinai diffusion [27, J8, _2£, 30] with disorder but no 
interaction corresponds to the limit c — > 0. 

Our focus is transport in this disordered, interacting 
particle, system. Since there is no overall bias in either 
direction, on average, the displacement of a particle with 
respect to its initial position, x, does not change with 
time, (x) = 0. We ask the most elementary question: 
how does the root-mean-square displacement, a, defined 



by cr^ 



evolve? We address this question via a 



scaling analysis of weakly disordered systems, e <C 1, and 
numerical simulations with general disorder strengths. 
Early Times. When disorder is weak, e <C 1, there is 
an initial period during which particles do not "feel" the 
disorder and hence move at random, = p_ = 1/2. In 
this early regime, disorder is irrelevant and the behav- 
ior is dominated by particle interactions. Without disor- 
der, the hard core repulsion causes a dramatic change 
in mobility: whereas an isolated particle moves diffu- 
sively, a ~ i^/^, the root-mean-square displacement of 
an interacting particle grows sub-diffusively with time 



laMEl, 22, 2, 



t 



1/4 



(1) 



Therefore, the movement of a particle is severely hin- 
dered by the presence of other particles. We illustrate 
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this remarkable collective behavior for extremely dense 
systems where there are large clusters of occupied 
sites that are separated by isolated vacancies. Particles 
move by exchanging their position with neighboring va- 
cancies. Furthermore, the sparse vacancies can be re- 
garded as non-interacting [20|. A particle that, up to 
time t, exchanges position with a total oi N — Af+ + 
vacancies of which iV+ were initially located to its right 
and were initially located to its left, has the dis- 
placement X = — First, since the vacancies are 
randomly distributed in the initial configuration, the ex- 
cess of vacancies in one direction follows from the law 
of large numbers, \N^ — N_ \ ~ N^/"^, and consequently, 
X ^ N^l"^. Second, vacancies that were initially located 
at a distance on the order of the diffusive length scale i^/^ 
from a particle may exchange position with it. There- 
fore, N ^ (1 — c) t^/^ and combining this scaling law 
with X ^ N^/^ yields ([T]). Although this scaling argu- 
ment applies to densely packed systems, the behavior ([1]) 
holds for arbitrary concentrations 0, H Hi- We 
also comment that this suppressed diffusion is a direct 
consequence of the hard core interactions. 
Intermediate Times. Eventually, the disorder be- 
comes relevant, and the biased hopping rates do affect 
the particle displacement. Although there is no global 
bias in the hopping rates, there certainly are local biases, 
as illustrated in figure[T]where sites with negative bias are 
in the majority. We expect that, at least at intermediate 
times scales, or equivalently, intermediate length scales. 
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these local biases lead to directed motion 

To quantify how such local biases affect particle mo- 
bility, we consider a particle that visits a distinct sites of 
which n_|_ have a positive bias and n_ have a negative bias 



with a 



Since the disorder is uncorrelated. 



the difference between the number of positive and nega- 
tive sites, A — — grows diffusively with the total 
number of visited sites, A ~ a^^^. The excess of sites bi- 
ased in one direction leads to a drift in this preferred di- 
rection with the small velocity v ~ e A/ct or v ^ ecr^^/^. 
Furthermore, the ballistic length scale x ^ vt gives an 
estimate for the displacement, x ~ (ei) cr~^/^. Since the 
displacement must be of the same order as the total num- 
ber of sites visited, a; ~ tr, we have 

-1/2 
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We thus arrive at our main result: the displacement be- 
comes super-diffusive because of the disorder, 

a^{etf\ (3) 

Of course, this length scale ultimately exceeds the sup- 
pressed diffusion length scale ([T]). We conclude that the 
inhomogeneous hopping rates generate a stochastic local 
velocity field, and as a result, there are local drifts that 
significantly enhance the mobility of the particles. 
Late Times. To understand the behavior at late times, 
we recall that the displacement of a non-interacting 



particle in a random disorder is logarithmically slow 
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At sufficiently large length scales, the random disorder 
generates a potential well that confines the particle. The 
depth of this potential well is the sum of all the biases in 
a given range, U{x) = X]i=i ^ ^^'^ there- 

fore, the depth of the well grows diffusively with distance, 
U ~ e^/x. This stochastic well constitutes a barrier that 
the particle must overcome, and since the time to escape 
out of this barrier grows exponentially with the depth of 
the well, t ~ exp(J7) ^ exp(ey^), the displacement is 
logarithmic as in ([4]). 

We argue that the slow mobility ^ also characterizes 
the asymptotic late time behavior of interacting particles 
in disorder. First, the confining potential well remains 
the same even when there are multiple particles. Second, 
the probability that a given particle escapes the well is 
exponentially small, and therefore, only mildly affected 
by the presence of other particles. We envision a sce- 
nario where particles are stuck in a local minimum of the 
potential and escape the barrier one at a time. Such an 
escape process is dominated by the same exponential es- 
cape time that characterizes an isolated, non-interacting 
particle. We conclude that at late times, interacting par- 
ticles in a random disorder also follow the logarithmic 
displacement law (j?]). Particle interactions become irrel- 
evant and the behavior is governed by disorder alone. 
The Three Time Regimes. By combining the early 
dl]), intermediate ([3]), and late @ time behaviors, we con- 
clude that the mobility of a given particle exhibits three 
distinct regimes of behavior as follows (see also figure [5]), 
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The time and length scales that characterize the crossover 
points can be obtained by matching the two correspond- 
ing behaviors. The transition from the early regime into 
the intermediate regime occurs at time t ~ e~8/5 ^^^^ 
length a ^ e~^/^, while the transition from the interme- 
diate domain into the late domain occurs at time t ^ e~^ 



34| and length a ^ e ^ , as shown in figure [51 



Let us recap the three regimes of behavior. At the 
early stages, particle interactions dominate over disor- 
der, and the motion of particles is sub-diffusive due to 
the hard core repulsion. In the intermediate regime, dis- 
order and interactions are both relevant. The particles 
stream following the stochastic local velocity field and 
the result is a strong, super-diffusive transport. At late 
times, disorder dominates and interactions become irrele- 
vant. Particles are trapped by a stochastic potential well 
and the displacement is logarithmically slow because the 
escape time is exponentially large. 
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t FIG. 3: The early and intermediate behaviors for weak dis- 

order, e = 10~^ (bullets) and no disorder, e = (squares). 
FIG. 2; The three regimes of behavior The displacement Shown is the displacement o versus time t (bullets), as well 
a is plotted versus time t using a double logarithmic scale. as a reference line with slope 2/3. 



As further support of the scaling behavior above, we 
can show that the stochastic potential well plays no role 
in the intermediate regime. Clearly, since the overall hop- 
ping rate equals one, the time scale characterizing the 
movement between neighboring sites is also of order one. 
The time to escape out of a well grows exponentially with 
the depth of the well, t ~ exp(C/), but this time scale be- 
comes appreciable only when the depth of the potential 
well is large, e\/x ^ 1, or equivalently, when the dis- 
placement becomes sufficiently large, x ^ e~^. Indeed, 
this length scale is realized only at the late time regime, 
as shown in figure [2l Therefore, trapping is negligible 
throughout the intermediate regime. 

Let us now consider the effect of disorder on a non- 
interacting particle. In the absence of disorder, e = 0, 
the particle displacement is unhindered and thus, purely 
diffusive, a ~ t^/'^. In weak disorder, e <C 1, an iso- 
lated particle undergoes ordinary diffusion at early times, 
but is later slowed down considerably according to ([3]). 
Hence, there are two distinct regimes of behavior when 
interactions are absent 
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We note that the crossover time scale matches the up- 
per time scale in ([5]). Thus, in the absence of particle 
interactions, disorder slows the particles down. 

Surprisingly, disorder has the opposite effect on an in- 
teracting particle system. Due to disorder, particles un- 
dergo fast, super-diffusive motion as in ([3]) over a sub- 
stantial time range and this effect becomes stronger as 
the disorder weakens because the crossover time and 
length scales are divergent. According to (O, the dis- 
placement in a given disorder eventually overtakes the 
displacement in a weaker disorder. This non-monotonic 
dependence on the disorder strength is another nontriv- 
ial consequence of the competition between disorder and 



interactions. We conclude that in the presence of particle 
interactions, disorder speeds the particles up. 
Numerical Simulations. We performed extensive 
Monte Carlo simulations to test the scaling predictions. 
The simulations are a straightforward implementation of 
the transport process. Initially, identical particles ran- 
domly occupy the sites of a one-dimensional lattice of 
size L with periodic boundary conditions, and the initial 
concentration equals c. Each lattice site has a bias in 
the positive or the negative direction as = 1/2 -I- e or 
p_ = 1/2 — e with equal probabilities. The dynamics are 
asynchronous. In an elementary step, a randomly chosen 
particle hops to the right with probability p+ or to the 
left with probability p_ , and this hop is successful only if 
the neighboring site is vacant. Subsequently, time is aug- 
mented by the inverse number of particles, t ^ t -{- 1/N. 
This elementary step is repeated indefinitely. We present 
results of simulations with L = 4 x 10^ and c — 1/2. 

We verified the super-diffusive behavior ^ using a 
weak disorder (figure [3]). Even though the super-diffusive 
behavior is an intermediate asymptotic, the duration of 
this regime grows rapidly as the disorder weakens. We 
performed a few additional tests: (i) we checked that the 
displacement cr is a function of the scaled time variable 
e t rather than t at intermediate times by using different 
disorders, (ii) we verified that the concentration does not 
play an important role using c = 1/4 and c = 3/4, and 
(iii) we used a different type of disorder with p+ drawn 
from a fiat distribution in the range [1/2 — e : 1/2 -I- e] 
and obtained qualitatively similar results. 

To test the behavior at late times, we also simulated 
a non-interacting particle system by ignoring the site oc- 
cupancy restriction. These simulations show that after 
an extremely long transient period, the displacements of 
interacting particles and non-interacting particles nearly 
match (figure thereby confirming that hard core in- 
teractions become irrelevant asymptotically, and that the 
behavior is governed by disorder alone. 
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ling transport properties because weak disorders result 
in strong mobilities. This effect can be tested experi- 
mentally in colloidal [l^l or biological 2^ 3| channels, 
where the slow transport ([T]) was realized. 
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FIG. 4: The late time behavior for e = 10~ . Shown is 
the displacement a versus time t for non-interacting particles 
(squares) and for interacting particles (bullets). 




FIG. 5: The behavior at moderate disorders. Shown is the 
displacement <j versus time t without disorder (e = 0, squares) 
and with moderate disorder values of e = 0.1 (bullets), 0.2 
(diamonds), 0.3 (down-triangle), and 0.4 (up-triangle). 



Our scaling analysis tacitly assumes that disorder is 
small. A comparison of the behaviors with moderate dis- 
orders and with no disorder shows that, irrespective of 
the disorder strength, mobility is always strengthened by 
disorder (figured]). Thus, mobility enhancement is a gen- 
eral effect that does not require weak disorder. 

Figure [5] also shows that the displacement in a homo- 
geneous system eventually catches up with the displace- 
ment in a strongly inhomogeneous system. Indeed, the 
sub-diffusive behavior ([1]) that characterizes a uniform 
system eventually exceeds the logarithmic displacement 
(|3|) in a disordered system. However, the crossover time 
t ^ is astronomical at weak disorders and in prac- 
tice, disorder always generates a stronger mobility. In- 
deed, the crossover time is very large even at moderate 
and strong disorders (figure [5|). 

In conclusion, we studied how disorder affects trans- 
port in an interacting particle system. We found that 
whereas disorder slows down non-interacting particles, 
disorder speeds up interacting particles. Therefore, there 
is an intricate interplay between interaction and disorder. 

Disorder provides an effective mechanism for control- 
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